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Abstract 

Three-dimensional (3D) turbulence has both energy and helicity as inviscid constants of 
motion. In contrast to two-dimensional (2D) turbulence, where a second inviscid invariant — 
the enstrophy — blocks the energy cascade to small scales, in 3D there is a joint cascade of 
both energy and helicity simultaneously to small scales. It has long been recognized that the 
crucial difference between 2D and 3D is that enstrophy is a nonnegative quantity whereas 
the helicity can have either sign. The basic cancellation mechanism which permits a joint 
cascade of energy and helicity is illuminated by means of the helical decomposition of the 
velocity into positively and negatively polarized waves. This decomposition is employed in 
the present study both theoretically and also in a numerical simulation of homogeneous and 
isotropic 3D turbulence. It is shown that the transfer of energy to small scales produces a 
tremendous growth of helicity separately in the -I- and — helical modes at high wavenumbers, 
diverging in the limit of infinite Reynolds number. However, because of a tendency to restore 



1 



reflection invariance at small scales, the net helicity from both modes remains finite in that 
limit. Since energy and helicity are not separately conserved in the + and — modes, there 
are four "flux-like" quantities for both invariants, which correspond to transfer either out of 
large-scales or into small scales and either to -I- helical or to — helical modes. The helicity 
fluxes out of large-scales in the separate + and — channels are not constant in wavcnumber 
up to the Kolmogorov dissipation wavenumber but only up to a smaller wavenumber kn , 
recently identified by Ditlevsen and Giuliani |^. However, contrary to their argument, 
the net helicity flux is shown to be constant all the way up to the Kolmogorov wavenumber: 
there is no shorter inertial-range for helicity cascade than for energy cascade. The transfer 
of energy and helicity between + and — modes, which permits the joint cascade, is shown 
to be due to two distinct physical processes, advection and vortex stretching. 
PACS numbers: 47.27.Ak,47.27.Gs 



1 Introduction 

Helicity is the signature of parity-breaking in incompressible fluid flows. For example, the 
most general two-point velocity correlation consistent with homogeneity and isotropy but not 
reflection symmetry is easily seen to have the Fourier transform 

/ *n +\ n +\\ fx kikj\ E{k,t) . k H{k,t) 

where E{k, t), H{k, t) are functions of the wavenumber magnitude k = |k| only. 5ij = 1 if i = j; 
0, otherwise. Siji = 1, if are cyclic; -1, if are anticyclic; 0, otherwise. It follows 

immediately from (1.1) by integration over k, both directly and after taking a cross-product 
with ik, that 

/•CO 1 POD 

/ dk E{k,t) = -{\v\'^) and / dk H{k,t) = {vu) , (1.2) 
Jo 2 Jo 

where u = Vxv is the vorticity field. Thus, E{k,t) is the standard energy spectrum, but 

H{k, t) is the spectrum of the helicity: 

H{t) = J dx v{x,t)-u{x,t). (1.3) 



It is only rather recently that it was realized that this quantity, like energy, is a quadratic 
invariant for the inviscid limit of the fluid equations, the incompressible Euler dynamics |3|, |^. 
Nonzero mean values of the helicity are now known to occur naturally in a wide variety of 
geophysical flows. For examples, tornadoes and updraft regions in rotating thunderstorms have 
helicity densities approaching 10 m/s^ and 0.1 m/s^, respectively [^]. For a review, see |^]. 

Since three-dimensional (3D) incompressible Navier-Stokes dynamics thus possesses two 
quadratic inviscid constants of motion, kinetic energy and helicity, it is natural to compare 
it with the situation for two-dimensional (2D) incompressible Navier-Stokes dynamics, which 
possesses also two quadratic invariants, the kinetic energy and the enstrophy (w^). As was 
shown by Kraichnan Q, the presence of the additional enstrophy invariant imposes a strong 
constraint on 2D turbulence, effectively blocking the energy cascade from large scales to small 
scales. As a result, there is in 2D a dual cascade, with energy cascading to large length-scales 
and with enstrophy cascading to small length-scales. The corresponding situation in 3D was 
first considered by Brissaud et al. who considered two possibilities for helicity cascades 
when the large-scales of the flow break reflection symmetry. The first possibility was a pure 
helicity cascade to small length-scales, in which, similar to the 2D case, there would be only 
a nonzero helicity flux 5 and no energy flux. In that case, Kolmogorov scaling arguments 
give energy and helicity spectra E{k) ~ 5'^/^k~'^/^, H{k) ~ 5'^/'^k~'^/^ . Note that when these 



spectra are substituted into (1^), the second parity- violating term is the same magnitude as 
the first term and never becomes negligible. The second possibility raised by ||8| was of a joint 
cascade in which there would be both a non-zero flux of energy e and also a non-zero flux of 
helicity 5 together to small length-scales. In that case, an assumption that the transfer rates 
are determined by e alone leads to spectra E{k) ~ e 

2/3^-5/3^ iJ(A:) ~ 5e~^/^k~^/'^. With these 
spectra, the second term in ( |1.1| ) becomes smaller than the first term at high wavenumber, by 
a factor l/Zc, and reflection symmetry is asymptotically restored. 

It was argued theoretically by Kraichnan in |Q that the second possibility is more plausible, 
based upon consideration of inviscid equilibria and the fact that a non-vanishing energy spec- 



trum E{k) does not imply a non- vanishing helicity spectrum H{k). On the contrary, in 2D, the 
enstrophy spectrum i^(k) = k'^E{k) and flux of energy to high-wavenumber always implies a 
corresponding enstrophy flux larger by the factor k'^ at wavenumber k. Numerical simulations 
also support the second possibility. In |1C] an EDQNM calculation found no evidence of a 



pure helicity cascade and instead demonstrated the existence of a joint cascade. Note that the 
assumption on transfer rates made by Q was built into the closure. A direct demonstration 
of the joint cascade picture was first given by [|ll| in a 128'^ direct numerical simulation of the 
Navier-Stokes equation with a large-scale helical forcing. A high Reynolds number was achieved 
in that simulation with hyperviscosity. We have repeated this simulation with a 512'^ resolution 
and with ordinary viscosity and verified the result [|l2|. 

In a recent pair of very interesting papers |2|, Ditlevsen and Giuliani have made a num- 
ber of striking new suggestions concerning the cascade state of energy and helicity. In par- 
ticular, they have argued for the existence of the joint cascade, but with a shorter inertial 
range for helicity than for energy. They have pointed out the existence of a new length-scale 
^ ~ {y I S)^!"^ e^l'^ which becomes, at high Reynolds number, much larger than the Kolmogorov 
scale r\ ~ In their picture, the joint cascade state with constant fiuxes e, b and spectra 

E{k) ~ £2/3^-5/3^ _ ^^-1/3^-5/3 ^^^^y ^^^^^ fo^ 

a range of wavenumbers l/L^fc^l/^. 
Thereafter, only the energy fiux remains constant ~ e and the spectrum of energy remains 
-E(A:) ~ e^l'^k~^l'^ but the inertial cascade of helicity is disrupted by viscous dynamics. In it 
has been verified that the proposed scenario holds for a popular shell model of turbulence, the 
GOY model, which has both energy-like and helicity-like inviscid invariants. For 3D Navier- 
Stokes the authors argue for their conclusions making important use of a decomposition of the 
velocity field into helical waves In their paper the authors call for a numerical test of 
their proposed scenario in 3D incompressible fiuids by a direct numerical simulation. 

These works of Ditlevsen and Giuliani are the immediate motivators of the present paper. 
We believe that the helical wave decomposition (HWD) employed by them is an important tool 
in the study of the joint energy-helicity cascade. It allows the contrast between 2D and 3D 



turbulence to be made most sharply. We follow jl], ^ in making use of the HWD, in even a 
more general form [14, 15 1. However, theoretically our conclusions are at odds with the main 



conclusions drawn by Ditlevsen and Giuliani in f^, In particular, we do not agree that the 
inertial range for helicity should be any shorter than for energy. We find that their argument 
neglects important cancellations that occur for the conserved quantities, as opposed to their 
helical components which are not separately conserved. In fact, the situation is considerably 
richer and more complex for 3D Navier-Stokes than for the simple GOY model. We find that 
there are quantities to which the considerations of ||l|, ^ apply, but others to which they do not. 
In addition to our new theoretical analysis, we also carry out here the numerical study by 3D 
DNS that was requested in ||l| . By means of it we have verified all of our important theoretical 
conclusions. 

The structure of this paper is as follows. In Section 2 we study the spectral fluxes of 
energy and helicity in the helical decomposition. We briefly review in subsection 2.1 the helical 
decomposition in Fourier space that we employ. In subsection 2.2 we introduce the helical 
decomposition of the energy and helicity spectra and their balance equations. In subsection 
2.3 we relate the transfers in the separate helical channels to the transfers of the conserved 
quantities. In subsection 2.4 we study the total transfer integrated over wavenumber and 
investigate, in particular, the signs and magnitudes of transfers between + and — helical modes. 
One important conclusion established there is that transfers will be in a direction so as to restore 
parity-symmetry at high wavenumbers. In subsection 2.5 we strengthen our exact analysis 
in the previous subsection by a Kolmogorov scaling analysis, which leads to a more refined 
picture. Here we develop our main predictions for the spectral fluxes in the joint cascade state. 
In subsection 2.6 we present the results on spectral quantities from a 512^ direct numerical 
simulation with helical forcing. Next, in Section 3 we turn to an examination of the joint 
cascade in physical space, using a filtering approach to distinguish large-scales and small-scales. 
In the first subsection 3.1 we present the helical decomposition in physical space. In subsection 
3.2 we apply this decomposition to dynamics of filtered quantities. Here we illuminate the 



important physical mechanisms contributing to transfer between + and — helical modes. These 
are crucial to give the cancellations which permit the joint cascade. In subsection 3.3 we consider 
the relation of filtered quantities to the spectral ones defined earlier. In subsection 3.4 we apply 
a Kolmogorov scaling analysis to predict the average behavior of the novel quantities which 
appear in the filtering approach. In subsection 3.5 we present our numerical results on the 
filtered quantities from the 3D DNS. Finally, in Section 4 we summarize our conclusions. 



2 Spectral Transfer 

2.1 Helical Wave Decomposition 

In [0, ^ use has been made of an expansion of incompressible velocity fields into helical waves 
|13| . We shall just remind the reader here briefly of this decomposition. If the fluid flow 
is contained in a periodic box, then the velocity fleld and vorticity fleld may be expanded 
into circularly-polarized or helical waves h±(k)e*'^*^, where h± are orthonormal and ikxh-t = 



±|k|h± [13|: 



v(x,0=5]5:a,(k,t)h,(k)e^i'- (2.1) 

k s=± 

and 

a;(x,t) = E E s\k\as{k,t)hs{k)e'''-^ (2.2) 

k s=± 

(Note that [|l3| took h*(k)'hs'(k) = 26ss', but we normalize to omit the factor 2.) We then 
define := o-th-i-e*''''' so that v = v^+v^. For the vorticity we instead use the convention 
that u:^ = X^k |k|a±h±e*^''' so that u = u;'^ - . Thus, = ±Vxv=^. 

Because the + and — helical modes belong to distinct eigenspaces of the curl operator, they 
are always orthogonal. It follows that the quadratic invariants of energy and helicity may be 
written as E = E^ + E~ and H = — , where 

E^{t) = IJ |v^(x,t)|2, H^{t) = J d^x a;±(x,t).v±(x,t) (2.3) 
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By the Parseval equality 

^^(*) = ^El«^(k>*)P, H^{t) = ^Ma^{Kt)f. (2.4) 
^ k k 

We see that in contrast to the net heUcity, which may have either sign, the partial helicities 

are always nonnegative (with our conventions). However, it is only the combinations 

E = £"•" + E~ and H = — H~ which are conserved by the inviscid dynamics, and not the 

separate terms. This has important consequences for the spectral dynamics in the it helical 

channels, which we now consider. 

2.2 The Helical Decomposition of Spectra and Transfer 

We may define nonnegative spectra for the helical modes via 

E^ik,t) = E hvtipM-pMk - IpI) (2.5) 
P ^ 

and 

HHk,t) = l{uHpM-p))d{k - |p|). (2.6) 
p ^ 

Note that^(A;,t) = E+{k,t) + E-{k,t) and H{k,t) = H+{k,t) - H' {k,t) give the usual energy 
and helicity spectra. 

It is easy to derive for these modal spectra, the following balance equations 

dtE'^{k,t) =T^{k,t) -2i^k'^E^{k,t) + F^{k,t) (2.7) 

and 

dtH^{k, t) = T^{k, t) - 2vk^H'^{k, t) + F^{k, t) (2.8) 

Here 

F^{k,t) = j2{it{p)M-pmk-\p\) (2.9) 
p 

with ^ = V ioi X = E and ^ = 2a; for X = H. This gives, ioi X = E and X = H, the input 
spectra of energy and helicity, respectively, into the positive and negative components of the 
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velocity field, by the force acting at large scales. The it transfer functions are defined similarly 
as in |11|, by 

T^{k,t) = J2N^{p,t)6{k-\p\) (2.10) 
p 

for X = E,H, where N^{p,t) = — {^f {p)[vjdjVi]{—p)) , and ^ has the same meaning as above. 
This formula will be used for numerical computation of T^{k, t). It is also a simple consequence 
of ( p.Sp and ( |2.10| ) that the standard input and transfer spectra are obtained as 

FE{k,t)=F+{k,t)+F+{k,t), TE{k,t)=T+{k,t)+T+{k,t), (2.11) 

and 

FH{k,t)=F+{k,t)-F+{k,t), TH{k,t)=T+{k,t)-T+{k,t). (2.12) 
An alternative set of formulas that are theoretically useful follow from 

E^{k,t)=Y: l{\a^{p,t)\')6{k-\p\) (2.13) 
p ^ 

and 

H^{k,t) = J2 |p|(|a±(p,t)|2)5(A; - |p|). (2.14) 
p 

where are (-v/2 times) the helical amplitudes defined by WalefFe. See |13], Eq.(5). From 
these we derive the very important relation that 

H^{k,t) = 2kE^{k,t). (2.15) 

The well-known inequality for the net spectra that H{k, t) < 2kE{k, t) ^ follows immediately. 
We see that it modes are separately maximal helicity modes. Eq.( p.l5 ) is the closest 3D 



analogue of the relation Q{k,t) = k'^E{k,t) which in 2D implies blocking of forward energy 
cascade by enstrophy conservation M. 
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2.3 Relation to Transfers of the Inviscid Invariants 



The ± parts of (|23),(|2]|) may be combined using (|2.12| ) to yield the standard spectral 



balance equations for energy and helicity: 

dtE{k, t) = TE{k, t) - 2vk'^E{k, t) + t) (2.16) 

and 

dtH{k, t) = Tnik, t) - 2vk^H{k, t) + FE{k, t) (2.17) 

However, there is a crucial difference between the transfer functions of the conserved quantities, 
Th,Te, on the one hand, and the transfer functions T^,T^, on the other hand. The integral 
relations 

dkTxik,t) = (2.18) 





for X = E,H are the statement of conservation of energy and of helicity, respectively. However, 
such relations are not satisfied by T^,T^ separately. 

The important relation ( p.l5|) allows us to obtain both the transfer functions and 
in terms of Te,Th- Indeed, as previously noted, 



E+{k,t)+E-{k,t) = E{k,t) (2.19) 
and from H+{k,t) -H~{k,t) = H{k,t) using (pl5| ) 



E+{k,t) - E-{k,t) = H{k,t)/{2k). (2.20) 

These can be solved to give 



,±n..^_U^n..^^H{k,t) 



and thus by ( ^.15 ) again 



E^{k,t) = -lE{k,t)±^^], (2.21) 



H^(k, t) = i {2kE{k, t) ± H{k, t)) . (2.22) 
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These equations are well-known pq] . However, if we then take the time derivative and use the 
balance relations for energy, we find moreover that 



1 



and that 



1 



TEik,t)± 



TH{k,t) 
2k 



T^{k, t) = - [2kTEik, t) ± THik, t)] , 



(2.23) 



(2.24) 



in terms of Te,Th- Of course, similar relations can also be developed for the forcing spectra 

^ E^^ H- 

We may define altogether four kinds of "flux-like" quantities 



nx'^(A;) = - Z'dpTKp) 
Jo 



and 



J k 



(2.25) 



(2.26) 



for each X = E, H. Thus, Ii^''"{k) represents the flow of X out of it modes for the wavenumbers 



< k and Iix^{k) represents the flow of X into the it modes for the wavenumbers > k. These 
quantities are generally not equal for < and >. In contrast, the net fluxes 11^ = H^''^ + 
n-'"", = n+''' - n^'"" for M=< or >, satisfy 



U<{k)=U>{k) 



(2.27) 



for X = E,H, as a consequence of ( 2.18| ). Therefore, no distinction need be made between the 
< and > quantities for the net fluxes. However, for the + and — modes separately the < and 
> fluxes may be different. We shall see that this is indeed the case. 



2.4 Total Balances in the Helical Decomposition 



A number of important qualitative results may be derived from the expressions ( 2.23 ), (|2.24 ) 
for the total balances in the it modes over all wavenumbers. If we integrate the relations 
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and ( |2.8| ) over wavenumber, we obtain 

dt 

and 



±R^ -D% + F^ (2.28) 



^^ = i?^-Z^| + F| (2.29) 

where D%{t) := 2v dk k^E^{k,t) and D^{t) := 2u dk k'^H^{k,t) are the total dissi- 
pation of energy and helicity, respectively, in the + and — components, and where F^{t) = 
/q°° dk F^{k, t) and F^{t) = /q°° dk F^{k, t) are the total forcing inputs of energy and hehcity, 
respectively, in the + and — components. For Re,Rh we obtain the expressions 

Re := ± ^dk T^{k,t) = dk (2.30) 



^0 



and 



roo roo 

Rh ■■= dkT^{k,t)= dk kTE{k,t) (2.31) 
Jo Jo 

by integrating (|]2|),(|2]24|) over wavenumber and using ( |2.18| ). 



As ( |2.28D , p.29[ ) show, the + and — components of energy and helicity are not separately 
conserved (in the absence of forcing and dissipation). In general, there is no reason that 
Re = oi Rh = 0. Instead, there will be a nonzero transfer between + and — components, 
Re for energy and Rh for helicity, which exactly balance, so that total energy and helicity are 
conserved. This has an interesting consequence. For total energy and helicity, the forcing input 
and the dissipation must balance, in the steady state: 

De = Fe = e, Dh = Fh = d. (2.32) 

However, for the + and — components separately, the input and dissipation do not balance, 
and in general there is a discrepancy: 

D%-F^ = ±Re, D^-F^ = Rh. (2.33) 



The signs of Re,Rh can be inferred from ( 2.30D ,( p.31| ). If the large-scales contain mostly 



positive helicity, then Tni^k) < at low wavenumbers and ( |2.30| ) implies Re < 0', conversely, if 
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the large-scales contain mostly negative helicity, then Tnik) > at low wavenumbers and (2.30) 
implies that Re > 0. Therefore, if we input for example mostly positive helicity at large scales, 
then < and > . This means that some of the energy carried by the + modes 
at large scales will be converted to energy of the — modes during transfer, and will then be 
accounted by the dissipation rather than by D^. We see that quite general considerations 
imply that there is a tendency toward equalization of energy in the + and — components due 
to the nonlinear transfer between them. Hence, reflection-symmetry which is broken at large 
scales will tend to be restored asymptotically at small scales. 



A similar argument using the formula (|2.31 ) for Rh implies that Rh > always, because 



TE{k) > at high wavenumbers in any case. Therefore, viscous dissipation of helicity in + and 
— components in general both exceed their forcing inputs: > F^. As the energy cascades 
to high wavenumber in each compoment, a corresponding helicity is carried in each component 
~ 2kE^{k) at wavenumber k, which in fact grows to exceed the input helicity in that component 
at large enough k. This is the same mechanism that, in 2D turbulence, prevents energy cascade 
to high wavenumbers due to enstrophy conservation. However, helicity has both signs, so that 
in 3D the growth in H~^{k) is compensated by an equal growth in H^{k) at high wavenumbers, 
with the net helicity conserved. The consequence is that both and will exceed their 
inputs. The equal growth of helicity in + and — modes at high wavenumber again implies an 
asymptotic restoration of parity-symmetry. 

2.5 Kolmogorov Scaling Theory 

A clearer picture is obtained from a Kolmogorov scaling analysis for steady-state helical tur- 
bulence in 3D. A similar analysis was attempted by Ditlevsen and Giuliani in but our 
assumptions and conclusions differ in some important respects from theirs. Ignoring for the 
moment any intermittency corrections, we assume that 

E{k) ~ CES^^^k-^/^ (2.34) 
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and 



H{k) ~ C7H(Ve / )fc 



1/3x^-5/3 



(2.35) 



both for <^ k <^ kE ^ e^l^v~'^l^. kE is the standard Kolmogorov wavenumber. The joint 
5/3 spectra for energy and hehcity were first proposed in fs], ^. The main assumption behind 
(|2.35| ) is that the transfer time for hehcity at wavenumber k is the same as for energy, TH{k) ~ 
TE{k) ~ e^^/^fc^^/'^ p. Equations ( 2.34|) ,( ^l35 ) along with ( ^.21 ) imply the corresponding 



energy spectra for components, already noted in ||l|, |T^ : 



E^{k) ~ ici,e2/3fc-5/3 



H^{k)^CEe^'''k-^'^ 





0] 


2k 




:ra 










2k \ 





(2.36) 



(2.37) 



for k^ ^ k <^ kE- However, Ditlevsen and Giuliani in arrive at a conclusion in con- 
tradiction to our starting assumptions ( |2.34| ), (|2.35 ). They define an intermediate wavenum- 
ber kn ~ e'^'^l'^ (b jv)^!"^ which is ^ kE at high Reynolds number. We refer to ku here- 
after as the "Ditlevsen wavenumber". It is concluded in Q that only for "the inertial range 
K\= our ki\ < k < kjj there is a coexisting cascade of energy and hehcity where helicity follows 
a 'linear cascade' with a H{k) ~ k~^^'^ spectrum. In the range kn < k < kE the dissipation of 
helicity dominates with a detailed balance between dissipation of positive and negative helici- 
ties and the right-left symmetry of the fiow is restored." We agree essentially with the latter 
statement, but we argue that the joint cascade of energy and helicity and the spectrum (2.35) 
hold up to the Kolmogorov wavenumber kE, not up just to the Ditlevsen wavenumber kn- 

We demonstrate the consistency of our assumptions. According to ( 2.36| ), the total energy 
separately in the it modes j/^^ dk E^{k) both remain finite as ^ 0, but by (|2.37] ) the helicity 
in the separate it modes increase as 
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/ dk H^{k) ~ e^/^v 
Jo 



3/4,-1/4 



(2.38) 
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for small u. Thus, both of these latter integrals will diverge in the limit of infinite Reynolds 
number. However, the total helicity dk H(k) will remain finite. This agrees with our 
conclusions in the previous subsection that helicity in the it modes separately will be produced 
as a consequence of energy transfer to higher wavenumbers in those modes. We may also 
consider total integrated dissipations. We see that in each of the it modes separately 

rkE 



D^ = 2u [ ^ dk k^E^ik) ~ i^e2/3^4/3 ^ ^2.39) 
Jo 



whereas 



= 2v dk k^H^ik) ~ ue^/^k'f oc e'/^i^~^/\ (2.40) 
Jo 

diverging as u ^ 0. This conclusion was already reached in However, we see with our 

assumption ( 2.35| ) that 

Dh = 2u dk k'^H{k) ~ i^{5/e^/^)k't{^ oc S, (2.41) 
Jo 

balancing the helicity input from the force, as it must. The important contribution to this 
integral comes from k ^ kE- On the contrary, Ditlevsen and Giuliani in |^ argued that "the 
integral will not be dominated by contributions from kE but contributions from /c//." This 
statement ignores precisely the cancellation which occurs between + and — modes in ( 2.41| ). 



However, the Ditlevsen wavenumber does play an important role in the behavior of the 
flux quantities n^'^(A;) and n^'^(/c) defined earlier. We show that both of these fluxes have 
constant values over a range of wavenumbers k^ <^ k <^ kn (at least), but that the constants are 
different. In fact, if we assume that the forcing is confined to a finite range of low wavenumbers, 
then integrating the relation (|2.8|) over the wavenumber interval [0, k] yields 



= n^'<(fc) + 2u r dp p^H^{p). (2.42) 
Jo 

Similarly, integrating ( |2.8| ) over the wavenumber interval [/c, 00] yields 

j-k 

D% = n|'>(A:) + 2v dp p^H^{p). (2.43) 
Jo 
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Now 2u j/^dp p^H^{p) ~ 1/^2/3^7/3 f^i^g^i ^j^is integral oc 5 precisely when k ~ kn as defined 
above. Since are generally of the same order of magnitude as 5 = — F^, 

n±'<(A;)^F|, (2.44) 

for the range ki <^ k <^ kn, smaller than the conventional inertial range ki <^ k <^ kE- On 
the other hand, 

n±'>(fc)«Z?± (2.45) 

over the conventional range, because the error term only becomes comparable to when k 
approaches kE- As we have seen, > F^ in general, with becoming arbitrarily large as 
Reynolds number increases. For k > kfj and k > kE, respectively, both the < and > flux-like 
quantities decrease below their "inertial-range" values, F^ and D^. In fact, 11^'^ (/c) decreases 
from to 0, while 11^'^ (/c) decreases from F^ to the negative value F^ — = —Rh- 
Ditlevsen and Guliani argued in that net helicity flux 

Unik) = n+'<(A;) - n^'<(A;) = n+'>(A;) - n^'>(fc) (2.46) 

would also have a constant value in only the same range, k^ <^ k <^ kfj, which is smaller than 
the conventional inertial range, ki k <^ kE- They illustrated this phenomemon in [|| for the 
GOY shell model. However, their argument ignores the cancellation that takes place between 
+ and — components for Navier-Stokes. In fact, for net helicity flux, 

Fh = Unik) + 2iy dp p^H{p) (2.47) 
Jo 

and 

Dh = IiH{h) + 2u dp p^H{p). (2.48) 

JO 

and the latter integral, 2v Jq dp p'^H{p), is negligible for all wavenumbers k in the range kL <C 
k <^ kE- Thus, Fh ~ Unik) ~ Dh for all k in the conventional inertial range, as expected. 
Exactly the same arguments as made above for helicity, when applied to energy, show that 

F| = n^'< (/c) + 2u t dp p^E^ (p) (2.49) 
Jo 
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and 



D| = n^'>(A;) + 1v dp p^E^ip). (2.50) 
Jo 

Thus, for the range <^ k <^ kE, 

n±'<(fe)«F|, (2.51) 



but 



U^^>{k)^D^. (2.52) 



The situation is quite similar to that for hehcity, except constancy of these flux-hke quantities 
holds over the whole extent of the conventional inertial range and the discrepancy — 
does not grow with the Reynolds number, but tends to a finite limit. 

2.6 Simulation Results 

To order to test the various theoretical predictions in subsection 2.1 ~ 2.5, we have carried 
out direct numerical simulation of incompressible isotropic turbulence by solving the forced 
Navier-Stokes equations using a pseudospectral calculation with 512^ resolution. The system 
was forced by holding the kinetic energy fixed in the first two shells for the wavenumber |k| < 2.5 
|17|. Nonnegative helicity was added by using a representation of the helicity spectrum in |18 



H{k,t)= 2k'-[R(k',t)xI(k',t)] (2.53) 

k'GS(A:) 

where S{k) is the wavenumber shell {k' : k-{Ak/2) < |k'| < k+{Ak/2)} and R(k, t) and I(k, t) 
denote the real and imaginary parts of the Fourier-transformed velocity v(k, t), respectively. 
At each time-step we rotated I(k, t) for all modes |k| < 2.5 so that its angle with R(k, t) was 
kept at 90° with positive handedness. This method was used by us earlier in |T^. It guarantees 
a positive helicity spectrum in the low wavenumbers and H{k, t) < 2kE(k, t) in general. The 
microscale Reynolds number achieved in our simulation was Rx = 220. The data analyzed were 
obtained after the system reached statistical equilibrium. The values of inputs of energy and 
helicity were e = 0.1427 and 5 = 0.4415. 
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In Fig. 1 are plotted the spectral fluxes of energy and helicity. 
Clearly, there is about a decade of inertial range where these fluxes have a constant mean 
value. It can be seen also that the intervals of constant flux for energy and helicity have about 
the same range in Fourier space, contrary to the prediction of ^ but in agreement with our 
modified theory. To really test this, we should increase the Reynolds number and show that the 
constant flux ranges for both invariants grow with the Reynolds number in the same manner. 
Computational limitations unfortunately prevent us from raising the Reynolds number beyond 
the present value. However, further evidence that the two constant flux ranges have the same 
extent will be given below for a shell model, where our arguments about cancellation also apply. 

Let us next consider the behavior in the separate it helical modes. Some of the relevant pa- 
rameter values in our simulations are given in Table 1. The values of the dissipations D^,D^ 
were obtained from the integrals ( 2.39D and ( |2.40| ), respectively. The transfers Re,Rh we 



calculated from the integrals (p.30|) and ( ^.311) , respectively. With our forcing technique, the 
inputs F^,F^ cannot be similarly obtained, because the forcing spectra are not easily com- 
putable. Instead, we have obtained those values via the relations (|2^ ), i.e. = D+ ^ Re 



and F^ = — Rh- We wish to examine the predictions made in section 2.5 for the fluxes 
n^'<,n^'> defined in (|2. 25^ ,(^3^) for ± polarization and X = E,H. We have predicted both 



the extents of the inertial ranges for these quantities and the plateau values in these ranges. 

In Fig. 2 (a) and (b) we plot n|'> and n|'<. 
We have predicted that all of these fluxes shall be constant over the same range ki <^ k <^ kE 
in the limit of high Reynolds number. It seems from the figure, however, that n^'*" has a 
somewhat shorter range than the other three fluxes. This does not contradict our predictions, 
since those were asymptotic statements at high Reynolds number. More precisely, the prediction 
is that all of the fluxes shall have ranges which grow at the same rate in that limit. This does 
not rule out their having ranges with lengths that differ by a constant factor. In fact, it is 
easy to see why H^'^ appears to have a shorter range in Figure 2. Note from the table above 
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Table 1: Parameters of the Flow 



Energy 


Helicity 


= \j.\Jil6 


L)jj = 7.30 


D]^ = 0.0685 


= 6.86 


Re = -0.05152 


Rh = 6.67 


F+ = 0.1239 


F+ = 0.63 


p- = 0.01698 


= 0.19 


e = 0.1427 


6 = 0.4415 



that is an order of magnitude smaller than in our simulation. The dissipation term in 
( 2.49| ) is negligible compared to the input value F^ for wavenumbers k <^ ^ [F^ / e^^^u)^^'^ , 
respectively, for it helical modes. More precisely, the wavenumber k^ can be defined as that 
for which the dissipation integral reaches a fixed fraction, say, 99%, of the input F^. Although 

grow at the same rate with Reynolds number, we see that, with our parameter values, fc^ 
should be only about 1/5 as big as k~^. This agrees very well with the results in Figure 2. 

Next let us consider scaling ranges for helicity. In Fig. 3 (a)and (b) we plot H^''* and H^'^. 
We have predicted that H^''* are constant from the injection wavenumber k^ up to the Kol- 
mogorov wavenumber kE, while H^'^ are constant only up to the Ditlevsen wavenumber kn 
(within constant factors). We can indeed see that the ranges for n|'^ are shorter. More quan- 
titatively, note that kfj/kE ~ {kh/kE)'^^'^ and since k^/kE ~ 0.1 in our simulation, we calculate 
that kfj/kE is about 1/3. This is quite consistent with the ranges observed in Figure 3. In 
contrast to the energy, both H^'^ and H^'^ have ranges shorter than the conventional inertial 
range, and they are expected to grow at a different rate than the Kolmogorov wavenumber. 
Furthermore, close examination of Figure 3 shows that the range for n^'*~ seems even a little 
shorter still than the range for n+'<. Ar guing just as for the energy, we see that the dissipation 
term in ( p. 42 ) is negligible compared with F^ for wavenumbers k <C k^ ~ (F^/e^/^z^)^/"^, re- 
spectively, for lb helical modes. Because F^ is about 3 times smaller than F^ in our simulation. 
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we predict that k^j should be only about 2/3 as big as k^. This is consistent with Figure 3. 

To verify our prediction for the plateau values of the fluxes, in Fig. 4(a) and (b) we plot 
n^'> and n^'<. 

In Fig. 4 (a), we can see that Hg'^ drop from values very close to in the inertial range to 
zero in the high wavenumber range, while we see in Fig. 4 (b) that H^'^ decrease from values 
close to the forcing inputs in the inertial range to values ^Re in the small scales. Likewise, 
in Fig. 5(a)(b) we plot n^'> and n^'<. 

There H^'^ drops from a value close to in the inertial range to zero at high wavenumber, 
and n^'^ from constant values at low wavenumbers to —Rh at high wavenumbers. Notice that 
these constants are somewhat smaller than the forcing inputs F^, especially in — modes, which 
is different from the prediction in ( 2.44| ) . The reason is that the inertial range in our simulation 



is not long enough so that we can not neglect the dissipation term in ( |2.42| ) . By calculating the 



dissipation explicitly, we have verified that the observations are consistent with (|2.42| ). 

It is clear that all of our theoretical predictions are consistent with the reported numerical 
simulation results, both qualitatively and quantitatively. Figure 1 demostrates the prediction 
that energy and helicity fluxes are constant over ranges of the same extent. Our results con- 
tradict the predictions of Ditlevsen and Giuliani that helicity flux shall be constant only up to 
the wavenumber kfj ~ e~'^l'^ (b jv)'^^'^ |||]. Since only a short inertial range can be achieved in 
our 3D simulations, it is worthwhile to give further evidence in a model with a longer inertial 
range. We therefore consider a helical shell model, which, like 3D Navier-Stokes has both + and 



— helical modes and energy and helicity as inviscid invariants 2C]. We consider a SABRA 



version of this model [21|, because the scaling properties are a little better than in the original 



GOY-type models. The model we study is, precisely, 



a^n'U„+2('"n+l) + ^^n-l^n+l(^n-l) 

ckn-2K-iK-2\ - ^klu-^ + /+ (2.54) 
aknul_^2iK+i)* + &fcn-i<+i(n^_i)* 
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-cA;„_2n+_iU+_2j - i^k^u:^ + /„ (2.55) 
for a set of modes and u~, n = 1,...,N. Here the wavenumber kn = X"'ko for some 



A > 1. This model may be termed SABRA3 in the classification scheme analogous to [20|. If 
a + 6 + c = 0, then this model for u = f = Q conserves both an "energy" E and "helicity" 
where 

N , TV 

E{t) = E 2^\<\' + Hit) = E - K\'] (2.56) 

n=l n=l 

and X = is the second root of the quadratic ax^ + 6x + c = 0, in addition to x = 1. We 



consider a standard choice of parameters a, b, c given in Table II of [20|, with A = 2 and a = 1. 
Thus, the second invariant is truly "helicity-like" . 

Unlike the GOY shell model considered by [|| , the SABRA3 model considered here can show 
the same cancellations between + and — modes that we have argued to occur in Navier-Sokes. 
If the argument of Q for Navier-Stokes were correct, it would apply also to this helical shell 
model. To check their prediction, we could vary either v or 6 and observe the changes in the 
lengths of the constant flux ranges. In fact, it is easiest to change 6. We performed two 22-shell 
simulations for SABRA3 with v = 10^"^, one with helicity input 6 = 0.0015 and another with 
helicity input 6 = 0.0003, 5 times smaller. If the prediction of Q were correct, then going 
from the first to the second, the inertial range for energy flux should be unchanged but the 
inertial range for helicity flux should be 2 times shorter. Fig. 6 shows the helicity flux in the 
two simulations with different helicity inputs. 

In the inset are the normalized energy and helicity fluxes. As we can see very convincingly from 
this plot, none of the constant flux ranges has its length changed by changing helicity input 6. 
This is in agreement with our prediction, but in disagreement with the prediction of |||]. 
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3 Filtering Approach 



3.1 Helical Decomposition in Physical Space 

In ||l|, H use has been made of an expansion of incompressible velocity fields into helical waves, 
for flow in a periodic box or infinite space ||l^. In fact, the helical decomposition is an intrinsic 
decomposition and can be made for any flow domain A. It corresponds just to a decomposition 
of the velocity field v(x, t) into eigenmodes of the non-dimensionalized curl operator S = 
(_A)-i/2vx [0,00 Acting on solenoidal vector fields (such as an incompressible velocity 
field), this operator satisfies = /. Hence, its eigenvalues are ±1. We can therefore write 
V = v+ + v^, with Sv^ = ibv^. In fact, formally = P^v, where = ^(1 it S) are the 
projectors onto it helical subspaces. Likewise, we can decompose the vorticity field as a; = 
cj^ —UJ^, for Scj^ = ibo;^. Notice our sign convention for defining . Thus, uj^ = ibVxv^. 
If the fluid flow is contained in a periodic box, then these decompositions agree with those 
obtained earlier by expansion into helical waves. 

Because the + and — helical modes belong to distinct eigenspaces of self-adjoint operator 
S, they are always orthogonal. It follows that the quadratic invariants of energy and helicity 
may be written as E = + E~ and H = — , where 

E^{t) = \Jd^^ |v=^(x,t)|2, H^{t) = J d^x c^±(x,t).v±(x,t) (3.1) 

Note in particular that 

H^{t) = (v±, (-A)i/2v±) > 0. (3.2) 

Thus, the partial helicities are always nonnegative (with our conventions). 

^Let us make explicit how this operator is defined. Suppose that e„(x) are the eigenfunctions of the minus 
Laplacian —A in the domain A. With appropriate boundary conditions aUowing integration by parts (e.g. 
Dirichlet), the corresponding eigenvalues \„ are nonnegative. Since the eigenfunctions form a complete set 
for any reasonable domain, we can expand any vector function v G as v(x) = ^^c„e„(x) for some 

vector coefficients c„. In that case, the operator is defined by (Ev)(x) = Ve„(x) XCn/x/XT. This gives a 
constructive definition of the operator that could be used in practical calculations. 
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3.2 The Helical Decomposition for Filtered Quantities 

As we have just shown, the hehcal decomposition can be made in any domain intrinsically and 
its applicability is not restricted to flow in a periodic box. We consider here the dynamics in 
a general domain, without use of Fourier analysis. Additional insight can, in fact, be obtained 
by considering the transfer dynamics of helical modes in physical space. However, to discuss 
transfer, we must distinguish different scales of motion. To resolve the dynamics simultaneously 
in space and in scale, we shall employ the same filtering approach that is used in large-eddy 
simulation [^2|. In this approach a low-pass filtered velocity v^(x, t) with scales < £ removed is 
introduced by a convolution vi = Gi * v with a smooth filtering function G£(x) = £~^G{'K/i). 
The filtered velocity obeys the equation 

9tv^(x,t) + V-[v^(x,t)v^(x,t) +TKx,t)] =fKx,t) - Vp^(x,t) + zy Av^(x,t) (3.3) 

Here ie,Pi are the filtered forcing and pressure, respectively, and Tg = (vv)^ — v^v^ is the tur- 
bulent stress, or spatial momentum transport induced by the eliminated small-scale turbulence. 
We may similarly consider a filtered vorticity a7^(x, t) and its corresponding equation (which is 
most simply obtained by taking the curl of ( |3.3| ) above). The total energy in the large-scales 
> i is represented in terms of filtered quantities by 

E{e,t) = ^J d^\vei^,t)\^. (3.4) 

and the total helicity in the large-scales > i hy 

H{i, t) = Jd^ v^x, t) ■ uJiix, t). (3.5) 

Note that these quantities represent cumulative values in the large-scales, not the analogues of 
Fourier spectra but instead spectra integrated in wavenumber from zero to 1 /i. 

We may assume the filtering operation to commute with S and (e.g. choose Gi = 
e^"^^). It follows that (v='=)^ = v^. Hence, we may introduce without any ambiguity the helical 
decompositions of large-scale velocity, V£ = v^ + v^, and of large-scale vorticity, = cjj' — uJJ . 
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The energy contained in it modes at length-scales > i is 

E^{i,t) = ^l d^\vfiK,t)\^ (3.6) 

so that E{i, t) = E^{i, t)+E~(£, t). Likewise, the helicity contained in it modes at length-scales 
> £ is 

H^{i,t) = J dx v±(x,t) •LZJ±(x,i). (3.7) 

so that H{£,t) = H~^{i,t) — H~{£,t). As we have defined these quantities, E'^{i,t) > and 
H^{£,t) > for both it helical modes. 

It is not difficult using the Navier-Stokes dynamics to obtain the evolution equations of 
large-scale energy and helicity. For example, for energy one easily calculates that 
dE^ 



dt 

where 



-(^, t) = t) - n|(£, t) ± RE{i, t) - D%{i, t) (3.5 



F|(^,t) = l^x f^(x,t).v±(x,t) (3.9) 
is the input of energy in the it modes at large length-scales > I, 

D^{e,t) = u J d:)^ |Vv^(x,t)|2 (3.10) 

is the dissipation of energy in the ± modes at length-scales > i, 

n|(^,t) = - fdx Vv^(x,t):T£(x,t) (3.11) 



is the net flux of energy from the it modes at length-scales > £ into the subgrid length-scales 
< i, and, finally, 

Re{£, t) = ±Jdx v7.(vrV)v± (3.12) 

is the transfer of energy entirely in the large length-scales between -|- and — components. We 
have already studied in ||l^ the statistics of the flux quantities Ii^{£,t). The formula ( 3.12 ) 



gives some further insight into the mechanism of transfer of energy between + and — modes. 
It is seen that a it-polarized helical wave, when advected by a mode of any other polarity, will 
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not stay pure it-polarized but will develop a component of the opposite sign polarity. This is 
the essential process contributing to RE{^-,t). 

Likewise, for helicity, it is not difficult to show that 



dt 

where 



, t) = F^ii, t) - n± (£, t) + Rnii, t) - D%ii, t) (3.13) 



F^ii,t) = 2 j d^ it{x,t)-Ljfix,t) (3.14) 
is the input of helicity in the it modes at large length-scales > i, 

D%{l,t) = 2u J dx Vv^(x,i):VaJ^(x,t) (3.15) 

is the dissipation of helicity in the it modes at length-scales > i, 

n|(^, = -2 / rfx VaJ±(x, ty.Teix, t) (3.16) 

is the net flux of helicity from the it modes at length-scales > i into the subgrid length-scales 
< i, and, finally, 

RH{i,t) = j [v7.(vrV)aJ± + v±.(vrV)aJ^ ± v±.(ZJrV)v^]. (3.17) 

Integration by parts shows that the latter expression is the same for both it signs. It represents 
a transfer of helicity entirely in the large length-scales between + and — components. We 
see that there are two mechanisms of such transfer. First, as for energy, advection of a it- 
polarized helical wave will produce a component of the opposite sign and this gives a transfer 
of helicity between -|- and — modes. In addition, a it-polarized helical wave which stretches 
a vortex of any polarity will generate vorticity of the opposite =F polarity. This mechanism of 
transfer of helicity between -|- and — modes is represented by the third term in (3.17). The 



vortex-stretching mechanism for transfer between + and — modes was already considered by 



Waleffe |13], who showed that it is the only mechanism that survives in the limit of large scale 



separation. 
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3.3 Relation Between Filtered and Spectral Quantities 

The quantities defined in the filtering approach above, when ensemble- averaged and with I set 
= 0, coincide with the corresponding total integrated spectral quantities defined in Section 2.3. 
Thus, from ( |3.8D ,( |3!T^ ) by taking the limit ^ ^ one rederives the equations ( p.28 ),( 2.29| ). Of 



course, in that limit the fluxes n^(£, t), n^(£, t) must vanish, for > 0. From this rederivation 
of (g38D,(|23g|), we obtain interesting identities for integrated transfers: 



oo 







dk ZMM = ±(vT.(v.V)v±), (3.18) 

Ak 



with both sides equal to Re^ and 

) 

dk kTE{k, t) = (v^-(v V)a;=^ + v^-(v V)cj^ ± v^-(a;-V)v^), (3.19) 



f 

Jo 



10 

with both sides equal to Rh- 

For ^ > 0, the equations (|3.8D , (|3.13| ) are essentially the same as the spectral balance equa- 



tions (2.7),( |2^ ) integrated in wavenumber from zero to ~ Clearly, it is the combined 



expression 

n^'<(£, t) := n|(£, t) ± Re{£, t) (3.20) 
which corresponds to the flux quantity Ii^'^{k,t) that we defined spectrally. Likewise it is 

n^'<(£,t) ■.= u%{e,t)-RH{i,t). (3.21) 

which corresponds to 11^'^ (A;, t). Just as for the spectrally defined quantities, 11^'^ (^,t) repre- 
sent the flow of X out of it modes for the length-scales > i, with X = E, H . There are two 
mechanisms that contribute to this flux. In the first, represented by Il^{i,t), X = E, H , the 
conserved quantity X leaves the it modes at length-scales > i and enters the subgrid modes with 
length-scales < i. In the second, represented by Rx{^,t), X = E, H, the conserved quantity 
X leaves the it modes at length-scales > £ hy being transferred to opposite polarity =p modes 
still at length-scales > £. The fluxes n^'^(^, t), n^'^(^, t) and the similar spectral quantities 
represent the net effect of both mechanisms. 
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3.4 Kolmogorov Scaling Theory 

The same Kolmogorov scaling arguments made in Section 2.4 in wavenumber space apply also 
to the filtered quantities. Hence, in the driven steady-state, the fluxes Il^'^{£,t),Il^'^{i,t) 
must have averages which are constant in the inertial-range: 

{U^'^'m^F^, (3.22) 

and 

(n±'<(£))«.F|, (3.23) 

as long as the mean dissipations (D^(^)), (D^(£)) can be neglected. This will hold for all 
L ^ i ^ rj = 1/kE in the case of n^'<(^) and for all L » ^ » e = ^/kn in the case of n^'<(£). 

However, the separate terms n^(£, i) and RE{i,t) should not be expected to have constant 
means. The average {RE{^,t)) will grow in magnitude through the inertial range, from near 
zero at ^ L to the value Re a^t very small i. For a very long inertial range, this should have 
value Re ^ ~ F^), as required to bring the energy in + and — modes into balance. 

The flux averages {11^(1)) will change from « ior i m L near the forcing scale to the value 

lb Re = -D^ for £ approaching the dissipation range. At high enough Reynolds number, 
the latter equals ^{F^ + F^) = | and, deep in the inertial range, the + and — modes should 
each carry about half of the total energy flux. Although (n|(^)) are not constant for £ moving 
through the inertial range, the net energy flux He = H^ + H^ satisfles 

{Ue{£)) « e (3.24) 

for all £ in the inertial range, where s = F^ + F^. Only the fractions carried by ± modes 
changes. 

Just as for energy, the separate terms n^(^, t) and Rh{£, t) are not expected to have constant 
means. In fact, for £ <^ L, one expects 

(H|(^)) ~ (3.25) 
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Such a result follows by the same reasoning that led Kraichnan to conclude for 2D that an 
energy flux at length-scale t would imply an enstrophy flux bigger by a factor ^"-^ at small- 
scales. This helicity flux continuously grows as I decreases through the inertial range, from a 
value ~ for i ^ L and eventually matches onto the dissipation ~ 2D^/r] oc e/ij at 
i rj. See ( p. 40 ). Of course, the growing term must cancel in the total flux — given by the 



difference H^'*^ = 11^ — Rh — which has constant mean F^, so that we have 

{RHm-''-^^^^-F^ (3.26) 

for I L m. the inertial range. This is nearly zero for £ ~ L, while taking i ^ r] implies 
that Rh ~ e/r/. This latter relation can be verified directly from the expression (|3l|). Using 
Tsik) = —-^IlEik) and integrating by parts gives 



oo 



Rh= dk UE{k) (3.27) 
Jo 

By Kolmogorov scaling analysis we expect that IlE{k) ~ e over a range of wavenumbers of 
length ~ kE- Thus, Rh ~ £ • kE- 

The divergence in the separate terms 11^,11^ cancels as well in the difference Hh = — 
n^;^, giving a constant net flux of helicity 

(UHie)) ~ 5 (3.28) 

for all £ ^ L in the inertial range, where S = F^ — F^. This means that the reflection symmetry- 
breaking at the large scales continues to be expressed in the statistics of the helicity flux 
throughout the inertial range. However, at sufficiently high Reynolds number and sufficiently 
deep in the inertial range, 

(n±(£))~^ + F±. (3.29) 
Therefore, parity-symmetry is restored to leading order. 
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3.5 Simulation Results 

We have checked the predictions in the preceding sections in our numerical simulation, as 
previously described. We employed a Gaussian filter, given in Fourier space as Gfc^(k) = 
g-(7r2k-k/24k2) ^-^^y^ _ Here, kc is the cutoff wavenumber and £ is the filter width. 

Our first result is the averages of the total filtered energy and helicity fluxes defined by 
IlE{i,t) = n^(^,t) + n^(^,t) and UH{e,t) = n^(£,t) - Ujj{£,t). These are plotted in Fig. 7, 
normalized by the total inputs e and S respectively. A short plateau between 30 < £/r) < 80 
indicates these filtered fluxes have a small inertial range. 

Let us now consider the fluxes in the ± helical channels, which we have argued not to be 
constant in the same range of scales. We begin with the energy. From the values reported 
earlier in Table 1, we expect that {Il'^{£))/e should decrease going to small scales through 
the inertial range, from /e ~ 0.90 at large scales to D~^/e ~ 0.53 at the beginning of the 
dissipation range. On the other hand, {Il'^{£))/e should increase through the inertial range, 
from /e 0.10 at large scales to D^/e ^ 0.47 at the start of the dissipation range. In the 
far dissipation range, both of the fluxes should decay to zero. The results in Fig. 8 agree well 
with these expectations. 

No theoretical prediction has been made for the rate of change of {n.^{£))/e through the inertial 
range, but they appear to be reasonably well fit by power laws shown in Fig. 8: 

{U+i£))/e^£^-'^ (n+(^))/£~r°-™. (3.30) 

Of course, this apparent power-law behavior can hold over a limited range only, else the limits 
D^/e at the beginning of the dissipation range could not be achieved. 

The increasing equahzation of (n^(£)) and (11^ (^)) at small length scales observed in Fig. 8 
is due to the transfer from — modes to + modes. It is expected that {Re{£))/£ will decrease 
monotonically from near zero at the beginning of the inertial range to a final value {Re)/£ ~ 
—0.37 in the far dissipation range (using the values in Table 1). This quantity is plotted in 
Fig. 9. 
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The expected behavior is weh-confirmed, including the asymptotic numerical value at very small 
scales. 

Now let us consider the helicity. It is predicted from ( 3.13| ) and ( |3.30| ) that {U'jj{i))/6 



should increase going to small scales through the inertial range, from ~ 1.43 at large 

scales to « 16.5 at the beginning of the dissipation range. In addition, {I[Jj(£))/6 should 

also increase through the inertial range, from ~ 0.43 at large scales to Djj/5 ~ 15.5. 

Proceeding into the far dissipation range, both of the fluxes should decay to zero. In Fig. 10, 
we show (n^(£))/5 and iJiJj{())/5 as functions of t/r]. 

Given the approximate power-law behavior observed for (n^(^))/e, we expect also approximate 
power-laws for lJl^{l))/5, because of the relation ( p. 25 ). More specifically, using ( 3.30| ), that 
relation predicts that 

{Yl+{i))/5^r'-'\ {Iij,{^))/5^r'\ (3.31) 

From the plots shown in Fig. 10, we see the numerical data are consistent with the above 
power law relations. In a very long inertial range it is expected (3.29) is valid for both of these 



quantities, which should scale as ~ The data obtained indicates that our simulation does 
not reach such an asymptotic limit. 

The average helicity transfer between it modes, (-R//(£))/5, is expected to increase from 
near zero to {Rh)/S ~ 15.1 from large scales to small scales based on Table 1. The numerical 
results presented in Fig. 11: 

As shown in Fig. 11, a power law {Ri{{i))/6 i^^'^ can be observed. Again, if a very long 
inertial range exists the exponent should be -1 not -1.5. From Fig. 11, we obtain {Rh)/6 ~ 11.22 
in the far dissipation range which is somewhat smaller than expected. This discrepancy is 
possibly attributed to the "smearing" effect of the filter. 

Unlike the above quantities with increasing and decreasing behaviors, the fluxes (11^'^ (^)) 
and (n^'*~(£)) are expected to have averages that are constant over the inertial range. Precisely, 



based on (lO^ ), (^1% and Table 1, (n^'<)/e should be 0.90 in the inertial range and 0.37 in 
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the far dissipation range, while (H^;' )/e should be 0.10 in the inertial range and -0.37 in the 
far dissipation range; similarly, (n^'^)/5 should be 1.43 in the inertial range and (n^'^)/5 
should be 0.43, while both go to —15.1 in the far dissipation range. As seen in Fig. 12, these 
expectations are well confirmed. 

Since the plateaux are short, the obtained values are somewhat smaller than expected. On the 
other hand, the values in the far dissipation range are quite close to those predicted. 

4 Conclusions 

A major theoretical goal of this work was to explain the mechanism which permits a joint for- 
ward cascade of energy and helicity in 3D turbulence. By contrast, a dual cascade of energy and 
enstrophy occurs in 2D turbulence, with a forward cascade of enstrophy and an inverse cascade 
of energy. It has long been understood that this difference arises from the sign-indefiniteness 
of helicity ||8|, ^, whereas enstrophy is nonnegative. We have argued that the contrast is seen 
most clearly in the helical decomposition of the 3D flow field, used previously in Just as has 
been argued in 2D for energy and enstrophy fluxes, in 3D a flux of energy to high wavenumber 
carries with it a growing flux of helicity in both -|- and — channels. The mechanism in 3D 
that allows the constant fluxes of energy and helicity to coexist is the near cancellation of the 
helicity flux between the -|- and — modes. The fact that these helical-channel contributions 
almost cancel is due to the asymptotic restoration of parity symmetry at high wavenumbers. In 
fact, we have shown by very general considerations that the nonlinear transfers of energy and 
helicity between + and — modes will be in such a direction as to restore reflection-invariance at 
small length-scales. In addition, we have made detailed, novel predictions about the constant 
flux ranges in the + and — channels, both their extents in wavenumber and the plateau values 
achieved. We have considered not only spectral fluxes, appropriate to a periodic box, but also 
fluxes in physical space for arbitrary space domains using a filtering technique. The predictions 
have been confirmed by a numerical simulation of forced helical turbulence. 
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The role of helicity in three-dimensional turbulence is, in our opinion, still somewhat mys- 
terious. In particular, it is still unclear how energy and helicity dynamics interact in detail. 
The role of helicity in geophysical flows has been considered Q — without being fully resolved — 
while its appearance and influence in engineering applications is still largely unexplored. We 
hope that this work will be a helpful step in the direction of better understanding the subtle 
manifestations of helicity in three-dimensional turbulence. 
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Figure Captions 

FIG. 1. Normalized energy and helicity flux in k-space. 
FIG. 2. (a)n^'<, n-'> and ,(b)n+'<, n+'> and Hg. 
FIG. 3. (a)n^'<, n^'> and UH,{h)U^'^, n^'> and Uh- 

FIG. 4. (a)n+'>, n~'> and (b)n+'<, n~'<. D+ = 0.0723, = 0.0685, i?E = -0.05152, 
= 0.1239, = 0.01698. 

FIG. 5. (a)n+'>, n^'> and (b) n+'<, n^'<. D+ = 7.30, = 6.86, Rh = 6.67, F+ = 0.63, 
= 0.19. 

FIG. 6. Normalized helicity flux in Helical SABRA3 model(22 shells). The inset shows the 
normalized energy (O) and helicity fluxes(«). 

FIG. 7. Normalized flltered energy and helicity flux. 

FIG. 8. Normalized filtered net energy fluxes in + and - modes, the solid lines are the power 
fitting. 

FIG. 9. Normalized filtered energy transfer between - and + modes. 

FIG. 10. Normalized filtered net helicity fiux in (a) + modes, (b) - modes, the solid lines are 
the power fitting. 

FIG. 11. Normalized filtered helicity transfer between - and + modes, the solid line is the 
power fitting. 

FIG. 12. (a) {U^'<{e,t))/e, {h){U%'<{e,t))/S. 
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